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Abstract 

We establish the optimal diversity -multiplexing (DM) tradeoff for coherent selective-fading multiple- 
access MIMO channels and provide corresponding code design criteria. As a byproduct, on the 
conceptual level, we find an interesting relation between the DM tradeoff framework and the notion 
of dominant error event regions, first introduced in the AWGN case by GaUager, IEEE Trans. IT, 1985. 
This relation allows us to accurately characterize the error mechanisms in MIMO fading multiple-access 
channels. In particular, we find that, for a given rate tuple, the maximum achievable diversity order 
is determined by a single outage event that dominates the total error probabUity exponentially in SNR. 
Finally, we examine the distributed space-time code construction proposed by Badr and Belfiore, Int. 
Zurich Seminar on Commun., 2008, using the code design criteria derived in this paper. 

I. Introduction 

The diversity -multiplexing (DM) tradeoff framework introduced by Zheng and Tse allows 

to efficiently characterize the information-theoretic performance limits of communication over 
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the STREP project No. IST-026905 MASCOT within the Sixth Framework Programme of the European Commission. Parts of 
this work were presented at the IEEE Int. Symp. Inf. Theory (ISIT), Toronto, ON, Canada, July 2008. 
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multiple-input multiple-output (MIMO) fading channels both in the point-to-point [1] and in the 
multiple-access (MA) case [2]. For coherent point-to-point flat-fading channels, DM tradeoff 
optimal code constructions have been reported in [3]-[6]. The optimal DM tradeoff in point-to- 
point selective-fading MIMO channels was characterized in [7]. In the MA case, the optimal 
DM tradeoff is known only for flat- fading channels [2]. A corresponding DM tradeoff optimal 
code construction was reported in [8]. 

Contributions. The aim of this paper is to characterize the DM tradeoff in selective-fading 
MIMO multiple-access channels (MACs) and to derive corresponding code design criteria. As 
a byproduct, on the conceptual level, we find an interesting relation between the DM tradeoff 
framework and the notion of dominant error event regions, first introduced in the case of 
additive white Gaussian noise (AWGN) MACs by Gallager [9]. This relation leads to an accurate 
characterization of the error mechanisms in MIMO fading MACs. Furthermore, we extend the 
techniques introduced in [7] for computing the DM tradeoff in point-to-point selective-fading 
channels to the MA case. Finally, we examine the distributed space-time code construction 
proposed in [10] using the code design criteria derived in this paper. 

Notation. Mt and Mr denote, respectively, the number of transmit antennas for each user and 
the number of receive antennas. The set of all users is W = {1, 2, ...,[/}, «S is a subset of W with 
S and |»S| denoting its complement in U and its cardinafity, respectively. The superscripts ^ and 
^ stand for transposition and conjugate transposition, respectively. A B and A B denote, 
respectively, the Kronecker and Hadamard products of the matrices A and B. If A has the 
columns {k—l, 2, . . . , m), vec(A) = [a{ af' ■ ■ ■ a^]^. ||a|| and || A||p denote, respectively, 
the Euclidean norm of the vector a and the Frobenius norm of the matrix A. For index sets 
Xi C {1, 2, . . . , n} and X2 C {1,2,..., m}, A(Xi, X2) stands for the (sub)matrix consisting of 
the rows of A indexed by Xi and the columns of A indexed by X2. The eigenvalues of the n x n 
Hermitian matrix A, sorted in ascending order, are denoted by Afc(A), k = 1, 2, . . . , n. The 
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Kronecker delta function is defined as 6n,m = 1 forn = m and zero otherwise. If X and Y are 
random variables (RVs), X r^Y denotes equivalence in distribution and Ex is the expectation 
operator with respect to (w.r.t.) the RV X. The random vector x ~ 0^(0, C) is zero-mean 
jointly proper Gaussian (JPG) with E{xx^} = C. f{x) and g{x) are said to be exponentially 
equal, denoted by f{x) = g{x), if lim^-^oo ^^^^ — l™a;-^oo ~bf^' Exponential inequality, 
indicated by > and <, is defined analogously. 

II. Channel and signal model 

We consider a selective-fading MAC where U users, with Mt transmit antennas each, com- 
municate with a single receiver with Mr antennas. The corresponding input-output relation is 
given by 



Yn = J^^^iiu,n^u,n + Zn: n = 0, 1, . . . , - 1, (1) 

u=l 

where the index n corresponds to a time, frequency, or time-frequency slot and SNR denotes 
the per-user signal-to-noise ratio at each receive antenna. The vectors y„, x„ „, and z„ denote, 
respectively, the Mr x 1 receive signal vector, the Mt x 1 transmit signal vector corresponding 
to the itth user, and the Mr x 1 zero-mean JPG noise vector satisfying E{z„z^} = 5„ „/ Im^, 
all for the nth slot. We assume that the receiver has perfect knowledge of all channels and the 
transmitters do not have channel state information (CSI) but know the channel law. 

We restrict our analysis to spatially uncorrected Rayleigh fading channels so that, for a given 
n, llu,n has i.i.d. CAA(0, 1) entries. The channels corresponding to different users are assumed 
to be statistically independent. We do, however, allow for correlation across n for a given u, and 
assume, for simplicity, that all scalar subchannels have the same correlation function so that 
E{H„,„(i, j) {Ilu',n'{h j))*} = RH(n, , for i = 1, 2, . . . , Mr and j = 1, 2, . . . , Mt. The 

covariance matrix Re is obtained from the channel's time-frequency correlation function [11]. 
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In the sequel, we let p = rank(Re). For any set S = {ui, . . . , u\s\}, we stack the corresponding 
users' channel matrices for a given slot index n according to 

H5j}2 [H^i-^^fi ■ ■ ■ Hui^ijJi]- (2) 

With this notation, it follows that 

E{vec(H5,„) (vec(H5,„/))^} = Re(^, n) I|5|MtMr- (3) 



III. Preliminaries 

Assuming that all users employ i.i.d. Gaussian codebook^ the set of achievable rate tuples 
i?2, • • • , Ru) for a given channel realization {iiu^n}n=o given by 
7^= I {RuR2,...,Ru):yS CU, 



fi(5)<lglogdet(l + ^H,„H|.„)} 



(4) 



where R{S) = J2ues ^ given rate tuple R2, . . . , Ru) ^ T^, we say that the channel 

is in outage w.r.t. this rate tuple. Denoting the corresponding outage event as O, we have 



P((9)=P U Os] (5) 
where the 5-outage event Os is defined as 

Os = < {lis,n}n=0 ■ 

> (6) 

-J^logdet I + — H,,.HjJ<i?(5) 

Our goal is to characterize ([5| as a function of the rate tuple {Ri, R2, . . . , Ru) in the high- 
SNR regime and to establish sufficient conditions on the users' codebooks to guarantee that the 

'a standard argument along the lines of that used to obtain [1, Eq. 9] shows that this assumption does not entail a loss of 
optimality in the high SNR regime, relevant to the DM tradeoff. 
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corresponding error probability is exponentially (in SNR) equal to F(0). To this end, we employ 
the DM tradeoff framework [1], which, in its MA version [2], will be briefly summarized next. 
In the DM tradeoff framework, the data rate of user u scales with SNR as i?u(SNR) = 
log SNR, where r„ denotes the multiplexing rate. Consequently, a sequence of codebooks 
Cr„(SNR), one for each SNR, is required. We say that this sequence of codebooks constitutes 
a family of codes Cr„ operating at multiplexing rate r„. The family C^^ is assumed to have 
block length A^. At any given SNR, Cr„(SNR) contains codewords = [x^ o ^u,i ■ ■ ■ ^u,n-i\ 
satisfying the per-user power constraint 

||X„||^<MTAr, VX„Ga„. (7) 

In the remainder of the paper, we will say "the power constraint (|7])" to mean that (|7]) has to be 
satisfied for u = 1,2, ... ,U. The overall family of codes is given by C^. = C^.^ x x ■ ■ ■ x C^^ , 
where r = {ri,r2, . . . , ru) denotes the multiplexing rate tuple|^ At a given SNR, the correspond- 
ing codebook Cr,(SNR) contains SNR^'^^^ codewords with r{U) = Y!i=i ^u- 
The DM tradeoff realized by Cr is characterized by the function 

diCr) = - lim ^^^^ 

^ SNR-.00 log SNR 

where Pe{Cr) is the total error probability obtained through maximum-likelihood (ML) decoding, 
that is, the probability for the receiver to decode at least one user in error. The optimal DM 
tradeoff curve d*(r) = sup,^^ d{Cr), where the supremum is taken over all possible families of 
codes satisfying the power constraint (|7]), quantifies the maximum achievable diversity order as 
a function of the multiplexing rate tuple r. Since the outage probability F{0) is a lower bound 
(exponentially in SNR) on the error probability of any coding scheme [2, Lemma 7], we have 

i.M<_ Ito ^ (8) 

^ ^ - SNR^oo log SNR ^ ^ 

^Throughout the paper, we consider multiplexing rate tuples lying within the boundaries determined by the ergodic capacity 
region. 
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where the outage event O, defined in ^ and is w.r.t. the rates i?,u(SNR) = log SNR, Vn. As 
an extension of the corresponding result for the flat-fading case [2], we shall show that ([8]) holds 
with equality also for selective-fading MACs. However, just like in the case of point-to-point 
channels, a direct characterization of the right-hand side (RHS) of ([8]) for the selective-fading 
case seems analytically intractable since one has to deal with the sum of correlated (recall that 
the H^j „ are correlated across n) terms in ([6]). In the next section, we show how the technique 
introduced in [7] for characterizing the DM tradeoff of point-to-point selective-fading MIMO 
channels can be extended to the MA case. 

IV. Computing the optimal DM tradeoff curve 
A. Lower bound on ¥{Os) 

First, we derive a lower bound on the individual terms ¥{Os)- We start by noting that for any 



set S ^U, Jensen's inequality provides the following upper bound: 

iV-l 

N 



^El°gdet(l + ^H,,.H-.) <logdet(l+l^^^^^^ (9) 



Us = < 



(10) 



where the "Jensen channel" [7] is defined as 

( 

[H5,o H5,i ■ ■ ■ H5,Ar_i], if Mr < |5|Mt, 

[H|o Hf,i ■ ■ ■ Hf,jv-i], if Mr > |5|Mt. 
Consequently, lis has dimension m(5) x A^M(5), where 

m(5) = min(|5|MT,MR) (11) 
M(5) = max(|5|MT, Mr). (12) 

In the following, we say that the event Js occurs if the Jensen channel Ti^ is in outage w.r.t. 
the rate r(5) log SNR, where r (5) = E«g5^«' i-^-' - (J^l^NR) < r(5)logSNR}. From 
we can conclude that, obviously, P(J5) < P(C5). 
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We next characterize the Jensen outage probability analytically. Recalling ([3]), we start by 
expressing Hs as Hs = ?iu,(R^/^ ® Im(5)), where R = Re, if Mr < \S\Mt, and R = R^, 
if Mr > |iS|Mt, and Tiu, is an i.i.d. CJ\f{0, 1) matrix with the same dimensions as Tis given by 



Tiw = < 



[H^,o H^,i ■ ■ ■ H^,jv-i], if Mr < \S\Mt, 

(13) 



[H^,oH^,i ■■■ a^-il ifMR>|5|MT. 
Here, iiw,n denotes i.i.d. CJ\f{0, 1) matrices of dimension Mr x |iS|Mt. Using ?iu,U ~ Tiw, 
for any unitary U, and A„(Re) = A„(Rj^) for all n, we get 

nsn^ ^u^{\® Im(5))^J (14) 

where A = diag{Ai(Re), AalRe), • • • , Ap(Re), 0, . . . , 0}. Setting Un, = : m(5)], [1 : 

pM(iS)]), it was shown in [7] that f{J's) is nothing but the outage probability of an effective 
MIMO channel with pM(S) transmit and m(iS) receive antennas and satisfies 

F{Js) =F(^\ogdet(l + SmUji!^)<r{S) logSNP) 

= SNR-'^^("('5)) (15) 

where we infer from the results in [1] that ds{r) is the piecewise Unear function connecting the 
points (r, ds{r)) for r = 0, 1, . . . , m(5), with 

ds{r) = (m(5) - r)(pM(5) - r). (16) 



Since, as akeady noted, F{Os) > P( it follows from ([B]) that 

P(C5) > SNR-^^^'^-^)) (17) 
which establishes the desired lower bound. 
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B. Error event analysis 

Following [2], [9], we decompose the total error probability into 2^ — 1 disjoint error events 
according to 

Pe{Cr) = J2 ^(^-5) (18) 

scu 

where the 5-error event £s corresponds to all the users in S being decoded in error and the 
remaining users being decoded correctly. More precisely, we have 

Ss = {{±u ^ X„, Vix e 5) A (X.„ = X,, e 5)} (19) 

where and X^ are, respectively, the transmitted and ML-decoded codewords corresponding 
to user u. We note that, in contrast to the outage events Os defined in ([6]), the error events Ss 
are disjoint. The following result establishes the DM tradeoff optimal code design criterion for 
a specific error event £3. 

Theorem 1: For every u E S, let Cr„ have block length > p|iS|Mt. Let the nonzercj^ 
eigenvalues of R;^0(^„g5 EfE„), where E„ = X„-X; andX„, X^ e C.JSNR), be given— 
in ascending order — at every SNR level by ?i„(SNR), n = 1,2, . . . , p\S\Mt. Furthermore, set 

m{5) 

A:lgT^(SNR)^ min J] ^^(^NR)- (20) 

{tLfu—J^u-J^uiueS ,_. 

X„,x;6C.„(SNR) "-^ 

If there exists an e > independent of SNR and r such that 

AmST^(SNR) > SNR-('^('5)-^), (21) 

then, under ML decoding, F{Ss) < SNR^'^^^''^'^^^ 

Proof: We start by deriving an upper bound on the average (w.r.t. the random chan- 



nel) pairwise error probability (PEP) of an iS-error event. From (19), we note that E„ 



[e«,o e„,i • • • eu,N-i], with e„,„ = x„,„ - x^„, is nonzero for m e 5 but E^ = for any 

^Here, we actually mean the eigenvalues that are not identically equal to zero for all SNR values. This fine point will be made 
clear in the proof. 

June 26, 2009 DRAFT 



u E S. Assuming, without loss of generality, that 5 = {1, . . . , \S\}, the probability of the ML 
decoder mistakenly deciding in favor of the codeword X' when X was actually transmitted can 
be upper-bounded in terms of X — X' = [Ei ■ ■ • E|5| ■ ■ ■ 0] as 

P(X X') 



SNR 

^ n=0 



(22) 



where 



Tl'(H5,„e„e«Hj„) 



Eh. 



with Us^n defined in Q and e„ = [e^^ „ ■ ■ ■ el^^^J^. Setting = [H^^o H^^i ■ • ■ H5,iv_i], 



we get from (22) 



P(X ^ X') 

< EH,|exp(^-^Tr (h^ diag{e„ef };;L"o Hf ^ 

= EH.{exp(-|^Tr (H^TT^H^)) } (23) 

where we have used = H^(R]^^^ I|5|Mt) with an Mr x A^|5|Mt matrix with i.i.d. 
CA/'(0, 1) entries and 

T = (R^/^ ® I|5|M,) dia.g{e^}^Z,' . (24) 

We note that T^T = (E«e5 E^E„), where we have rank(E„g5 E^E„) < |5|Mt be- 
cause Eu has dimension Mt x and > |5|Mt by assumption. Recalling that rank(Re) = p 
and using the property rank(A B) < rank(A) rank(B), it follows that rank(T^T) < 
p|(S|Mt, which is to say that T^T has at most p|iS|Mt eigenvalues that are not identically 
equal to zero for all SNRs. We stress, however, that these eigenvalues may decay to zero as a 
function of SNR. Next, using the fact that for any matrix A the nonzero eigenvalues of AA^ 
equal the nonzero eigenvalues of A^A, the assumption (made in the statement of the theorem) 
that R]^ (Xlue^ E^E„) has p|iS|Mt eigenvalues that are not identically equal to zero for all 
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SNRs implies that so does YT^ . The remainder of the proof proceeds along the lines of the 
proof of [12, Th. 1]. In particular, we split and subsequently bound the 5-error probability as 



+ F{Js)F{£s\Js) 

<i 

<nJs)+^{Ss\Js). (25) 



As detailed in the proof for the point-to-point case given in [12], the code design criterion (21 ) 



yields the following upper bound on the second term in ( 25 1 



- 7j^^ ]■ (26) 
In contrast to the Jensen outage probabihty which satisfies F{Js) = SNR"''^^''^'^^^ the RHS of 



(26l decays exponentially in SNR. Hence, upon inserting (26l into (25 ), we get ^{8^) < ¥{J's), 



and can therefore conclude that F{£s) < SNR^'^^^'^'^^^ 



In summary, for every £s, (21 1 constitutes a sufficient condition on {Cr„ : u E S} for F{£s) 
to be exponentially upper-bounded by F{J's). This condition is nothing but the DM tradeoff 
optimal code design criterion for a point-to-point channel with |iS|Mt transmit antennas and Mr 
receive antennas presented in [12]. In order to satisfy this condition, the users' codebooks have to 
be designed jointly. We stress, however, that this does not require cooperation among users at the 
time of communication. We are now ready to establish the optimal DM tradeoff for the selective- 
fading MAC and provide corresponding design criteria on the overall family of codes Cr. 

C. Optimal code design 
We start by noting that ^ implies F{0) > F{Os) for any S C U, which combined with 



(17) gives rise to 2^ — 1 lower bounds on F{0). For a given multiplexing rate tuple r, the 
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tightest lower bound (exponentially in SNR) corresponds to the set S that yields the smallest 
SNR exponent ds{r{S)). More precisely, the tightest lower bound is 

P(C>) > SNR-'^^^WS*)) (-27) 
with the dominant outage event given by Os*, where 

S* = aigmm ds{r{S)) (28) 
scu 

is the dominant outage set. Next, we show that, for any multiplexing rate tuple, the total error 



probability Pe{Cr) can be made exponentially equal to the RHS of (27 ) by appropriate design of 



the users' codebooks. As a direct consequence thereof, using ([8j), (27), and Pe{Cr) > f{0) [2, 
Lemma 7], we then obtain that ds*{r(S*)) constitutes the optimal DM tradeoff of the selective- 
fading MIMO MAC. Before presenting this result, let us define the function rs{d) as the inverse 
of ds{r), i.e., d = ds(^rg{d)^ and r = ^•^((/^(r)). We note that rg{d) is a decreasing function 
of d and ds{r) is a decreasing function of r. 

Theorem 2: The optimal DM tradeoff of the selective-fading MIMO MAC in ([!]) is given by 
d*{r) = ds^riS")), that is 

d*{r) = {m{S*) - r{S*)){pM{S*) - r{S*)). (29) 



Moreover, if the overall family of codes Cr satisfies ( pT[ ) for the dominant outage set S* and, 
for every S ^ S*, there exists e > such that 

^ST^(SNR) > SNR" (^5-^) (30) 

where 

0<ls<rsids*ir{S''))) (31) 

then 

d{Cr) = d*{r). (32) 
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Proof: Using (fTSl), we write 



(33) 



We bound the terms in the sum on the RHS of ( [33] ) separately. By assumption, Cr satisfies (21 ) 
for S* and, hence, it follows from Theorem [T] and ([15]) that 



(34) 



Next, we consider the terms F{£s) for S ^ S* and use (25 ) to write 



n£s)<njs)+^{£s\Js) 



(35) 



where (35 1 is obtained by the same reasoning as used in the proof of Theorem [T] with the users' 



codebooks {Cr^ : u E S} satisfying (30l instead of (21 1. Inserting (34) and p5\ into (33 ) yields 



= SNR^'^'S*^''^'^*)) 



(36) 
(37) 



where (37) follows from the fact that (31 ) implies ds{^s) > ds*{r{S*)), for all S ^ S*, and 



consequently, the dominant outage event dominates the upper bound on the total error probability. 



With Pe{Cr) > P(0) [2, Lemma 7], combining (27) and (37) yields 



PeiCr) = P(0) = SNR-'^5*("('5*)). 



(38) 



Since, by definition, d{Cr) < d*{r), using (^, we can finally conclude from (38) that 



d{Cr) = d\r) = dsMSn)- 



(39) 



As a consequence of Theorem [2] the optimal DM tradeoff is determined by the tradeoff 
curve ds*{r{S*)), which is simply the SNR exponent of the Jensen outage probability P(j75*) 
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corresponding to the dominant outage set. By virtue of ( 15 1, (38 1, and the fact that the relations 
F{Os) < P(C) and F{Js) < F{Os) hold for every S and, a fortiori, for the dominant outage 
set S*, we get 

nos*)<no) = njs* ) < nos* ) (40) 

which is to say that 

F{Os*) = SNR-'^5*("('^*)\ (41) 

Hence, as in the point-to-point case [12], the Jensen upper bound on mutual information yields 
a lower bound on the outage probability which is exponentially tight (in SNR). 

In order to achieve DM tradeoff optimal performance, the families of codes {Cr^,u E hi} 



are required to satisfy (21 ) for the dominant outage set S* and, in addition, the probability 



P(£^5) corresponding to the sets S ^ should decay at least as fast as ¥{Os*) = F{Js*), a 



requirement that is guaranteed when (30) is satisfied for every S ^ S*. Note that this code 



design criterion is less stringent than requiring all the terms F{Ss) to satisfy condition (21 1, as 
originally proposed in [13, Th. 2]. We conclude by pointing out that the code design criterion in 
Theorem |2] was shown to be necessary and sufficient for DM tradeoff optimality in Rayleigh 
flat- fading MACs in [14]. We stress, however, that there exist codes — at least in the two-user 
flat- fading case — that satisfy ( [2T] ) in Theorem [1] for all 5 CU as we will show in Section [V] 



D. Dominant outage event regions 

The following example illustrates the application of Theorem|2]to the two-user case, and reveals 
the existence of multiplexing rate regions dominated by different outage events. Remarkably, 
although the error mechanism at play here (outage) is different from the one in [9], the dominant 
outage event regions we obtain have a striking resemblance to the dominant error event regions 
found in [9]. 
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Example: We assume Mx = 3, Mr = 4, and rank(Re) = p = 2. For U = 2, the 2^ — 1 = 3 
possible outage events are denoted by Oi (user 1 is in outage), O2 (user 2 is in outage) and O3 
(the channel obtained by concatenating both users' channels into an equivalent point-to-point 
channel is in outage). The SNR exponents of the corresponding outage probabihties are obtained 



from ( [161 ) 

c?n(r„) = (3-r„)(8-r„), u = 1,2, (42) 
+ rs) = (4 - (n + rs)) (l2 - (n + rs)) . (43) 



Based on (42 1 and (43 1, we can now explicitly determine the dominant outage event for every 
multiplexing rate tuple r = (ri, r2). In Fig.[T] we plot the rate regions dominated by the different 
outage events. Note that the boundaries ri < 3, r2 < 3, and ri + r2 < 4 are determined by 
the ergodic capacity region. In the rate region dominated by Oi, we have c/i(ri) < d2{r2) and 
di{ri) < (i3(ri + r2), implying that the SNR exponent of the total error probability equals di(ri), 
i.e., the SNR exponent that would be obtained in a point-to-point selective-fading MIMO channel 
with Mt = 3, Mr = 4, and p = 2. The same reasoning applies to the rate region dominated 
by O2 and, hence, we can conclude that, in the sense of the DM tradeoff, the performance in 
regions Oi and O2 is not affected by the presence of the respective other user. In contrast, in the 
area dominated by O^, we have d3{ri + r2) < du{ru), u = 1,2, which is to say that multiuser 
interference does have an impact on the DM tradeoff and reduces the diversity order that would 
be obtained if only one user were present. 

Fig. |2] shows the dominant outage event regions for the same system parameters as above 
but with one additional receive antenna, i.e., Mr = 5. We observe that not only larger sum 
multiplexing rates are achievable, i.e., ri + r2 < 5, but also that the area where dominates the 
total error probability, and hence where multiuser interference reduces the achievable diversity 
order, is significantly smaller relative to the area dominated by the single user outage events Oi 
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Fig. 1. Dominant outage event regions for a two-user MA MIMO channel with Mt = 3, Mr = 4, and p = 2. 



and 02- This effect can be attributed to the fact that increasing Mr yields more spatial degrees of 
freedom at the receiver and, consequently, alleviates the task of resolving multiuser interference. 

E. Multiplexing rate region at a given diversity level 

The dominant outage event determines the maximum achievable diversity order as a function 
of the multiplexing rate tuple r. Conversely, one may also be interested in finding the region 
lZ{d) of achievable multiplexing rates at a minimum diversity order d G [0, pMtMr] associated 
with the total error probability. This is accomplished by designing an overall family of codes 
that is DM tradeoff optimal and satisfies 

ds{r{S))>d, ySCU (44) 
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Fig. 2. Dominant outage event regions for a two-user MA MIMO channel with Mt = 3, Mr = 5, and p = 2. 



which upon application oirs{-) to both sides is found to be equivalent to 

r{S)<rs{d), ySCU. 

We just proved the following extension of [2, Th. 2] to selective-fading MA MIMO channels. 

Corollary 1: Consider an overall family of codes Cr that achieves the optimal DM tradeoff 
in the sense of Theorem |2j Then, the region of multiplexing rates for which the total error 
probability decays with SNR exponent at least equal to d is characterized by 

7^(rf) = jr : r{S) < rs{d),WS C w| (45) 

where rs{d) is the inverse function of ds(r). 

To illustrate the concept of a multiplexing rate region [2], consider the two-user case with 
Mt = 3, Mr = 4, and p = 2. Fig.[3]shows the multiplexing rate regions TZ(d) corresponding to 
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several diversity order levels, i.e., d E {0, 2, 4, 8, 16}. The region 7^(0) is the pentagon described 
by the constraints ri < 3, r2 < 3, and ri + r2 < min(2MT, Mr) = 4. Higher diversity order 
can be achieved at the expense of tighter constraints on the achievable multiplexing rates ri 
and r2. For instance, for a diversity order requirement of c? > 8, the achievable multiplexing 
rate region is given by the pentagon OABCD. Increasing the minimum required diversity order 
results in multiplexing rate regions that shrink towards the origin. Note that to realize a diversity 
order requirement of ci > 16, the allowed multiplexing rate region is a square; in this case, 
performance (in the sense of the DM tradeoff) is not affected by the presence of a second user. 
Intuitively, the required diversity order is so high that users can only communicate at very small 
multiplexing rates and multiuser interference does not dominate the total error probability. 

V. Analysis of a code construction for the two-user flat-fading case 

In this section, we study the algebraic code construction proposed recently in [10] for flat- 
fading MACs with two single-antenna users and an arbitrary number of antennas at the receiver. 
We examine whether this code satisfies the code design criteria of Theorem |2] and focus on the 
case of a two-antenna receiver, for simplicity. 

We start by briefly reviewing the code construction described in [10] for a system with Mt = 1, 
Mr = 2,U = 2, N = 2, and p = 1 (i.e., flat fading). For each user u, let Au denote a QAM 
constellation with 2^"'^^'^'^) points carved from Z[z] = {k + il : k,l E Z}, where i = and 
<(SNR) = (r„ - e) logSNR for some e > 0, i.e., 

f _2K(SNR)/2 2^"(SNR)/2 ^ 

Au = Uk + U): <k,l< ,k,lEZ\. (46) 

The proposed code spans two slots so that the vector of information symbols corresponding to 
user u is given by s„ = [s^^i s^ ^], where s„^i, s„_2 ^ Au- The vector Su is then encoded using 
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Fig. 3. Multiplexing rate regions as a function of the diversity order d € {0, 2, 4, 8, 16} corresponding to the total error 
probability (Mt = 3, Mr = 4, and p = 2). 



the unitary transformation matrix U underlying the Golden Code [5] according to 



x; = u s„ 



a(x„ 



withU 



V5 



(47) 



where </? = denotes the Golden number with corresponding conjugate (p — a — 
1 + i — i(p and a — 1 + i — i(p. By construction, a;„ belongs to the quadratic extension Q{i, y/b) 
over Q{i) — {k + il : k,l E Q}. Here, a denotes the generator of the Galois group of Q(i, v^) 
given by 

a: Q{i,V5) ^ Q(i,V5) 
a + by/5 i— > a — by/5. 



(48) 
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Moreover, one of the users, say user 2, multiplies the symbol transmitted in the first slot by a 
constant 7 G Q{i), resulting in the transmit codeword 



X 



(49) 



Xi cr(xi) 
7x2 a{x2) 

Depending on the choice of the parameter 7, the codeword difference matrices arising from 
this construction have a nonzero determinant. For completeness, we shall next provide a proof 
of this statement which was originally made in [10]. 



Theorem 3: For any 7 7^ ±1 and any two X, X' according to (49 1, it holds that det(A) 7^ 0, 
where A = X - X'. 

Proof: Proceeding along the same lines as [10], we start by proving that the determinant 



corresponding to any codeword X in (49 1 is nonzero for any 7 7^ ±1, and hence, by the linearity 
of the mapping a over Q(i, v^), the determinant of any codeword difference matrix is also 
nonzero. Note that 



det(X) = Xia{x2) — 7X2cr(xi^ 



X — 7cr(x) 



(50) 



where the last step follows from setting x = Xia{x2), noting that a{a{x)) = x for any x G 
Q{i, V5), and using the property a{x ■ y) = a{x) ■ a{y) for every x,y E Q{i, V5). Hence, 
det(X) is zero if and only if 7 satisfies 7 = x/a{x). In this case, recalling that 7 G Q(^), 
we must have x G Q(i), or x G \/5Q(i) = {v^(A; + il) : k,l E Q}. These constraints yield, 
respectively, 7 = x/a{x) = 1 and 7 = x/a{x) = — 1, from which we can infer that det(X) = 
<;=^ 7 = ±1. Hence, any 7 G Q(2)\{±1} guarantees det(X) 7^ for xi,X2 E Q(i, V^). ■ 
We are now ready to examine whether this construction satisfies the code design criteria for 
DM tradeoff optimality given in Theorem |2j For simplicity, we assume 7 = Hn the following. 
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We start by considering the cases 5 = {1} and S = {2}. Assume that Au is chosen according 



to ([46|). By (47 1 and the fact that U is unitary, we obtain 



max = max ||Sn||^ 



2HJ,(SNR) 



(51) 



for M = 1, 2. In order to satisfy the power constraint (|7]), we scale the transmit vector corre- 
sponding to user u as 

/2i?4(SNR)\ -1/2 ^ 

x,„ = I z I x„ (52) 



so that, using (51 1, we get 



max Xt, 

X„ GCr„(SNR) 



2«;.(SNR) 



max ||Xt, 



2. 



(53) 



For user 2, we note that (53) remains valid after multiplying the first entry of X2 by 7 = i. 
Next, we note that in the flat-fading case Rj^ (e^e„) = e^e„, where = x„ — x^ for 
x„,x^ G Cr.„(SNR) and m = 1,2. Considering user 1, i.e., S = {!}, we have \S\ = 1 and 



m{S) = 1 so that the quantity defined in (20) is simply the smallest squared norm of the first 



row in (57 ) and satisfies 



Aj(SNR) = min 

XI ,x'^ eCri(SNR) 



I ' l|2 

|Xi - x^ll 



2i-^i(SNR) 

xi : XI = SiU^ ; x'^ : i[ = s'^U'^ 
Si,2, «i,i> «i,2 6 -^1 



|~ -/ ||2 




(54) 



(55) 
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where (54) follows from (52), and (55 ) is a consequence of = Us^ and the unitarity of U. 



From (46 ), we note that d^i^ = 1, i.e., the minimum distance in Ai is independent of SNR, and 



invoking _R'^(SNR) = (ri — e) log SNR, we can conclude from (55 ) that 



Ai(SNR) = SNR-("^~^). 
For user 2, a similar argumenj^ shows that Aj^(SNR) = SNR^*^*^^""^ and, hence, the code satisfies 



the criteria arising from (21 ) for 5 = {1} and S = {2}. 



Next, we consider the case S = {1,2}. The overall transmit codeword is now given by 

2-i?^{SNR)/2 2-«2(SNR)/2 

and satisfies the power constraint (|7]), i.e.. 



X = ^ 



(56) 



max IIX 

XgC,.(SNR) 



\l = max Tr (XX^) 

XeCr(SNR) ^ ^ 



II II 2 II II 2 

max Xi + max X2 

xieCri(SNR) xaSCraCSNR) 



From (56) and the linearity of the mapping a over Q(i, y^), the codeword difference matrix 



E = V2 



(57) 



is obtained as 

2-i?'i(SNR)/2 2--f^i(SNR)/2 ^^g^^ 

2-R^(SNR)/2 2-^2(SNR)/2 ^^g^^ 

where e„ = x„ — x'^ and hence e« G Q(i, v^), u = 1,2. Recall that in the flat-fading 
case Rj^ (E^E) = E^E. Next, note that \S\ = 2 and m{S) = 2 so that Ai(SNR) = 



minE |det(E) 1^. From (57 ), we readily get 



min |det(E)|' = 2i-(^'i(SNR)+^4(SNR)) 

E=X-X' 
X,X' e C^(SNR) 



mm 

A=X-X' 

X=/(si,S2),X'=/(si,s^) 



|det(A) 



(58) 



The multiplication of the first component of X2 by 7 = i does not affect the Euclidean norm. 
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where we have used the notation X = /(si, S2) to express the fact that X is obtained from Si 



and S2 using (47) and (49). We recall that det(A) 7^ for A arising from any combination 
of vectors s„, (m = 1, 2) with entries in Z[z]. Therefore, for every SNR, there must exist an 
cu(SNR) > such that 

min |det(A)|^ = Lu(SNR) (59) 

A=X-X' 

X=/{si,S2),X'=/(s'i,s^) 



which, upon inserting into (58) and using _R^(SNR) = (r^ — e) log SNR {u = 1, 2), yields 



A2(SNR) = min ^ |det(E)|^ 

X,X' e Cr(SNR) 



= SNR-("^+"^-2^)cu(SNR). (60) 

It follows from ( [59] ) — by inspection — that a; (S N R) is a nonincreasing function of S N R. Unfortu- 
nately, Theorem[3]does not allow us to conclude that (SNR) is bounded away from zero, in which 



case we could conclude from ( 60 ) and (21) that the code is DMT-optimal for all multiplexing rate 
tuples. Therefore, characterizing the decay of (SNR) as a function of SNR is key to proving or 
disproving the DM tradeoff optimality of the code construction. Unfortunately, we have not been 
able to determine how co'(SNR) decays with SNR|^ Characterizing this decay rate seems very dif- 
ficult and is likely to require advanced algebraic concepts. We can, however, distinguish between 
three different possibilities. If u;(SNR) decays exponentially with SNR, the criteria for DM 
tradeoff optimality provided in this paper are not met. In the case of a subpolynomial decay, i.e., 

'°ftr=o 

SNR-00 log SNR 

we would get A2(SNR) = SNR^*^'^^^'^^^^-' and, hence, such a decay would be sufficient to 



guarantee that (60) satisfies the code design criterion (21 ) for 5 = {1, 2} and any tuple (ri, r2 



^We would like to use this chance to point out that despite the claim we made in [13], we do not have a proof establishing 
the DM tradeoff optimality of the code construction in [10]. 
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in the multiplexing rate region. Finally, we consider the case of Ci;(SNR) exhibiting polynomial 
decay, assuming that ci;(SNR) = SNR^^, 5 > 0. In this case, it would follow from (60) that 



A^(SNR) = SNR-("i+"^+^-2^\ 



The quantity A|(SNR) would then decay faster than required by (21 1. In other words, the code 
construction would not be DM tradeoff optimal in the sense of Theorem [2] when the dominant 
outage set is 5* = {1,2}. However, when the dominant outage set is either S* = {1} or 
S* = {2}, the relaxed (compared to the code design criteria proposed in [13]) code design 



criteria provided in (30) would still be met for any multiplexing rate tuple (ri, satisfying 

ri + r2 + 5 < rs{ds*{r{S*))). 

We conclude this section by noting that a DM tradeoff optimal code construction for flat-fading 
MACs was reported in [8]. Specifically, it is shown in [8] that lattice-based space -time codes 
achieve the optimal DM tradeoff with lattice decoding. As a consequence of the code design 
criterion in Theorem [2] being necessary and sufficient for DM tradeoff optimality in Rayleigh 
flat-fading MACs [14], the code construction reported in [8] necessarily satisfies these design 
criteria. The systematic construction of DM tradeoff optimal codes for selective-fading MA 
MIMO channels seems, however, largely unexplored. 

VI. Conclusion 

We characterized the optimum DM tradeoff for selective-fading MA MIMO channels and 
studied corresponding code design criteria. Our results show that, for a prescribed multiplexing 
rate tuple, the optimal DM tradeoff is determined by the dominant outage event. The systematic 
design of DM tradeoff optimal codes for the (selective-fading) MIMO MAC remains an important 
open problem. 
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